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Abstract 

In  the  solid  oxide  fuel  cell  (SOFC)  system,  the  internal  reforming  of  raw  fuel  will  act  as  an  efficient  cooling  system.  To  realize  this 
cooling  system,  a  special  design  of  the  internal  reformer  is  required  to  avoid  the  inhomogeneous  temperature  distribution  caused  by  the 
strong  endothermic  reforming  reaction  at  the  entrance  of  the  internal  reformer.  For  this  purpose,  a  tubular  internal  reformer  with  adjusted 
catalyst  density  can  be  inserted  into  the  tubular  SOFC  stack.  By  arranging  this,  the  raw  fuel  flows  along  the  axis  of  the  internal  reformer  to 
be  moderately  reformed  and  returns  at  the  end  of  the  internal  reformer  as  a  sufficiently  reformed  fuel. 

In  this  paper,  the  output  characteristics  of  this  configuration  are  simulated  using  mathematical  models,  in  which  one-dimensional 
temperature  and  molar  distributions  are  computed  along  the  flow  direction.  By  properly  mounting  the  catalyst  density  in  the  internal 
reformer,  the  temperature  distribution  of  the  cell  stack  becomes  moderate,  and  the  power  generation  efficiency  and  the  exhaust  gas 
temperature  are  higher.  Effects  of  other  operating  conditions  such  as  fuel  recirculation,  fuel  inlet  temperature,  air  recirculation  and  air  inlet 
temperature  are  also  examined  under  the  condition  where  the  maximum  temperature  of  the  stack  is  kept  at  1300  K  by  adjusting  the  air  flow 
rate.  Under  this  condition,  these  operating  conditions  exert  a  considerable  effect  on  the  exhaust  temperature  but  have  a  slight  effect  on  the 
efficiency.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Tubular  solid  oxide  fuel  cells  (SOFCs)  have  been  pointed 
out  to  have  advantages  with  respect  to  mechanical  stresses 
developed  by  seals  and  thermal  stresses  caused  by  the 
temperature  distribution  along  the  cell  axis,  in  comparison 
with  planar  SOFCs  [1],  However,  leveling  of  the  tempera¬ 
ture  distribution  is  still  required  to  avoid  cell  deterioration 
due  to  the  higher  local  temperature.  The  higher  temperature 
may  lead  to  electrode  sintering  and  chemical  reaction 
between  the  electrode  and  the  electrolyte.  Therefore,  the 
temperature  distribution  has  to  be  restricted  in  a  narrow 
range. 

The  reforming  of  raw  fuel  is  an  endothermic  reaction. 
Therefore,  the  internal  reforming  can  be  used  to  cool  the  cell 
stack,  only  if  some  design  is  contrived  to  avoid  the  inho¬ 
mogeneous  temperature  distribution  caused  by  the  strong 
reforming  reaction  at  the  inlet  of  the  internal  reformer. 
To  fill  this  requirement,  the  internal  reformer  should  be 
separated  from  the  fuel  electrode  to  adjust  the  catalyst 
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density.  For  this  purpose,  a  tubular  internal  reformer  with 
a  thin  or  graded  catalyst  density  can  be  inserted  into  the 
tubular  cell.  The  raw  fuel  will  flow  then  along  the  axis  of 
the  reformer  to  be  reformed  moderately  and  will  change 
its  flow  direction  as  a  sufficiently  reformed  fuel  at  the  end 
of  the  internal  reformer. 

There  is  also  another  advantage  of  the  internal  reforming 
avoiding  the  heat  loss  that  cannot  be  avoided  in  the  external 
reforming. 

In  this  paper,  the  schematic  structure  of  the  tubular  SOFC 
equipped  with  an  internal  reformer  is  presented,  and  the 
effects  of  the  internal  reforming  and  other  operating  con¬ 
ditions  are  discussed  by  numerical  simulations. 


2.  SOFC  system  equipped  with  internal  reformer 

A  schematic  structure  of  the  SOFC  system  having  both  a 
pre-reformer  and  an  internal  reformer  is  shown  in  Fig.  1.  In 
this  figure,  one  cell  stack  with  an  internal  reformer  and  two 
air  supply  tubes  are  shown.  In  a  real  system,  many  stacks, 
internal  reformers  and  air  tubes  have  to  be  arranged.  The 
enlarged  drawing  of  the  stack  will  be  shown  later  in  Fig.  3. 
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Nomenclature 

Ac  activity  of  carbon 

Ag  cross-section  of  gas  CV  concerned  (cm2) 

Are f  reforming  area  per  length  (cnr/cm) 

As  cross-sectional  area  of  solid  CV  concerned 

(cm2) 

Asks  equivalent  value  considering  the  effect  of 

multilayers 

c  total  molar  concentration  (mol/cm3) 

Cp.  heat  capacity  of  gas  species  j  (J/mol) 

de  equivalent  diameter  (cm) 

Djj  diffusivity  between  gas  species  i  and  j  in 

porous  media  (cm2/s) 

Dy  diffusivity  between  gas  species  i  and  j  in 

free  flow  (cm2/s) 

E  electromotive  force  of  relevant  position  of 

cell  at  T  (V) 

E(i )  electromotive  force  of  CV(0  concerned 

E0(T)  standard  electromotive  force  at  T  (V) 

F  Faraday’s  constant  (96485  C/mol) 

AGact  activation  energy  (J/mol) 

/ig  heat  transfer  coefficient  of  mixed  gas 

(W/cm2  K) 

Hj{T)  enthalpy  of  gas  species  j  at  T  (J/mol) 
Hy(T)  heat  of  formation  of  gas  species  j  at  T  (J/mol) 

A//£H4  enthalpy  change  in  combustion  of  CH4  at 
298  K  in  HHV  (J/mol) 

A//j3oo  enthalpy  change  of  reaction  at  1300  K 
(kJ/mol) 

4(0  current  flowing  in  air  electrode  of  C V(0 

(A) 

4(0  current  flowing  in  electrolyte  of  CV(0  (A) 

4(0  current  flowing  in  fuel  electrode  of  CV(0 

(A) 

4  load  current  (A) 

J  current  density  in  electrolyte  (A/cm2) 

4  exchange  current  density  (A/cm2) 

kg  thermal  conductivity  of  gas  (W/cm  K) 

ks  thermal  conductivity  of  solid  (W/cm  K) 

Kc  equilibrium  constant  of  the  Bourdouard 

reaction 

Kox  equilibrium  constant  of  oxidation  of 

hydrogen 

Ks  equilibrium  constant  of  water  gas  shift 

reaction 

Mj  mole  flow  rate  of  gas  species  j  (mol) 

My  mole  flow  rate  of  air  gas  species  j  (mol) 

()/Wch4  mole  of  reformed  methane  in  CV  concerned 
(mol/s) 

SMco2  mole  of  C02  generated  by  water  gas  shift 
reaction  in  CV  concerned  (mol/s) 

/VceM  number  of  cells  in  stack 

Nj  mole  flux  of  gas  species  j  (mol/cm2  s) 

Pj  partial  pressure  of  gas  species  j  (atm) 


P,  B  partial  pressure  of  gas  species  j  at  boundary 

between  electrolyte  and  fuel  electrode  (atm) 
Po2,b  partial  pressure  of  oxygen  at  boundary 

between  electrolyte  and  fuel  electrode  (atm) 
^o2  partial  pressure  of  oxygen  in  air  (atm) 

^o2ib  partial  pressure  of  oxygen  at  boundary 

between  electrolyte  and  air  electrode  (atm) 
Gceii.f  heat  migration  to  side  wall  with  fuel  flow 
due  to  cell  reaction  (W) 

Gceii.ox  heat  migration  to  side  wall  with  oxygen 
flow  due  to  cell  reaction  (W) 

2conv(0  enthalpy  migration  due  to  gas  flow  to  CV(0 
(W) 

Gheat  heat  generated  in  CV  concerned  (W) 

Gjoui  heat  generated  by  Joule  heating  (W) 

Grad  net  heat  input  by  radiation  (W) 

Greact  enthalpy  migration  to  side  wall  due  to 

convection  (W) 

Greform  heat  generated  by  reforming  of  methane  (W) 
G shift  heat  generated  by  water  gas  shift  reaction  (W) 

Gcond  heat  input  due  to  heat  conduction  from 

neighboring  gas  CVs  (W) 

'fZjQj  summation  of  quantity  Q  in  relevant  gas 
flow 

mstQj  summation  of  quantity  Qj  in  exhaust  gas 
xpinputairgr  summati0n  of  quantity  Qj  in  input  air 
r  radial  distance  in  cylindrical  coordinate  (cm) 

rsub  radius  of  substrate  (cm) 

R  universal  gas  constant  (8.31441  J/(mol  K)) 

Rtl  branch  resistance  of  air  electrode  (Q) 

Rcih  rate  of  reforming  of  methane  (mol/(g  Ni  s)) 

Re  branch  resistance  of  electrolyte  (Q) 

Recf  recirculation  rate  of  fuel 

Rf  branch  resistance  of  fuel  electrode  (Q) 

/?int  branch  resistance  of  interconnecting  mate¬ 

rial  (£1) 

Rt  reforming  rate  of  methane 

Ss  side  wall  area  of  CV  concerned  (cm2) 

Ssl  side  wall  area  of  inside  CV  (cm2) 

Ss2  side  wall  area  of  outside  CV  (cm2) 

T,  T0  temperature  (K) 

4  temperature  of  air  (K) 

4 o  air  inlet  temperature  to  system  (K) 

4in  air  inlet  temperature  to  cell  (K) 

4eli(0  temperature  of  cell  CV(i)  (K) 

4eii,max  maximum  temperature  of  cell  stack  (K) 
4(0  temperature  of  fuel  CV(0  (K) 

4in  fuel  inlet  temperature  to  cell  (K) 

4o  fuel  inlet  temperature  to  system  (K) 

4(0  temperature  of  gas  CV(0  (K) 

4nax  maximum  allowable  temperature  (K) 

4x(0  temperature  of  oxidant  CV(0  (K) 

4m  exhaust  gas  temperature  from  SOFC  system 

(K) 

4(0  temperature  of  solid  C V(0  (K) 
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rsx  temperature  of  opposite  solid  CV  (K) 

Uf  fuel  utilization  rate 

t/ox  oxygen  utilization 

Vact  activation  polarization  (V) 

Kct,f  activation  polarization  of  fuel  electrode  (V) 

Kct,ox  activation  polarization  of  air  electrode  (V) 
yceii  average  value  of  cell  voltage  (V) 

Vcons  concentration  polarization 

yp  polarization  voltage  (V) 

Vp(0  polarization  voltage  of  CV(i)  (V) 

Vt  load  voltage  of  cell  stack  (V) 

x  flow  direction  coordinate  (cm) 

Sx  length  of  CV  concerned  (cm) 

Xj  mole  fraction  of  gas  species  j 

Xjn  mole  fraction  of  gas  species  j  at  boundary 

between  electrolyte  and  fuel  electrode 
*o2  mole  fraction  of  oxygen  in  oxidant 

*o2,b  mole  fraction  of  oxygen  at  boundary 

between  electrolyte  and  air  electrode 
Wcat  catalyst  Ni  weight  density  (g/cm2) 

Greek  letters 

a  transfer  coefficient 

[>  proportional  constant  (1/cm2  Q) 

e  porosity 

k  reaction  order  of  oxygen 

>)dc  thermal  efficiency  of  SOFC  in  HHV 

a0  interface  conductivity  (1/cm2  Q) 

Subscripts  and  superscripts 

B  boundary  between  electrode  and  electrolyte 

f  fuel 

g  gas  (fuel,  fuel  (reformed  fuel),  air,  oxidant) 

ox  oxidant 

s  solid  (tube,  cell,  internal  reformer) 


The  residual  spent  gases  are  burned  in  a  combustor  to 
obtain  the  higher  temperature  of  the  exhaust  gas  out  of  the 
SOFC  system.  The  combustion  gas  heats  the  pre-reformer 
and  the  supplied  air  sequentially. 

If  the  catalyst,  Ni,  is  loaded  densely  and  uniformly  on 
the  internal  surface  of  the  internal  reformer  in  the  manner 
used  for  the  fuel  electrode,  the  reforming  reaction  may  be 
completed  at  the  entrance,  and  this  region  may  be  cooled 
strongly  to  distort  the  whole  temperature  distribution.  To 
avoid  this,  the  density  of  the  catalyst  has  to  be  thin  in  order 
to  complete  the  reforming  reaction  at  the  end  of  the 
internal  reformer.  If  a  more  complicated  design  can  be 
allowed,  the  density  has  to  be  thinner  at  the  inlet  region 
and  thicker  at  the  outlet  region.  By  this  grading  of  the 
catalyst  density,  the  temperature  distribution  may  become 
flatter. 


The  raw  fuel,  methane,  comes  into  the  pre-reformer  at  the 
left  end,  then  goes  into  the  internal  reformer.  The  reformed 
fuel  changes  its  direction  and  flows  back  through  the  annulus 
formed  by  the  internal  reformer  and  the  cell  stack.  The  air 
comes  into  the  air  supply  tube  at  the  right  end  after  being 
preheated  by  the  exhaust  gas  through  the  heat  exchanger. 
The  air  also  changes  its  direction  at  the  left  end  and  flows 
back  as  an  oxidant.  Some  parts  of  the  spent  gases  are 
recirculated  using  ejectors  to  preheat  the  inlet  gases. 

The  recirculated  spent  fuel  also  supplies  the  steam 
required  to  reform  methane. 

The  cell  stack  is  a  power  generator;  on  the  other  hand  it  is 
a  heater.  In  the  SOFC  system,  the  cooling  of  the  cell  is 
carried  out  by  the  reactants  at  the  lower  temperature  and  the 
reforming  reaction.  The  flow  rates  of  the  reactants,  the  gas 
recirculation  rates,  the  inlet  temperatures  of  the  reactants 
and  the  rate  of  reforming  are  interdependent  under  the 
restricting  condition  to  maintain  the  maximum  temperature 
of  the  cell  stack  at  the  preset  value. 


3.  Definition  of  thermal  efficiency  of  SOFC  system 


The  thermal  efficiency  of  the  methane-fuelled  SOFC 
(abbreviated  “efficiency”  here),  rjdc  is  defined  as 

VtIo  _  keen  ,, 

^  AW-A//«hJ  1.1534^’ 

because 


where  Vt  is  the  dc  output  voltage  of  cell  stack,  I0  the  load 
current,  MCH4  the  mass  flow  of  input  methane,  A//®H4  the 
enthalpy  change  in  the  combustion  of  methane  at  298  K, 
Vceii  the  average  cell  voltage,  Uf  the  fuel  utilization  rate  and 
/Vcen  the  number  of  cells  in  the  stack. 

As  can  be  seen  in  Eq.  (1),  Uf  must  be  higher  to  obtain 
higher  efficiency;  however,  if  Uf  is  too  high,  the  fuel  will  be 
depleted  at  the  end  of  the  cell  and  the  cell  voltage  will  drop 
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substantially.  Therefore,  Uf  is  chose  to  be  85%  in  our  system 
referring  to  the  reported  experimental  results  [2 — 4]. 


4.  Numerical  models  for  simulation 


For  the  simulation  of  the  SOFC  system  with  the  inter¬ 
nal  reformer,  we  need  numerical  models  of  reforming  the 
cell  and  heat  migration.  These  models  will  be  described 
sequentially. 

In  gas  flow  modeling,  plug  flow  is  assumed  in  the  main 
flow  region,  and  flow  is  assumed  in  the  porous  media  such  as 
the  electrodes  and  the  substrate. 

A  set  of  the  cell  stack  equipped  with  an  internal  reformer 
and  an  air  injection  tube  will  be  considered  from  now  on, 
assuming  periodicity  in  the  F-axis  direction.  Namely,  the 
effect  of  the  heat  loss  from  the  vessel  is  ignored.  The  X-axis 
originates  from  the  inlet  of  the  internal  reformer  as  indicated 
in  Fig.  1. 

4.1.  Reforming 

In  the  reforming  of  methane  with  steam,  the  dominant 
reactions  are  the  following  two  reactions: 

CH4  +  H20  =  3H2  +  CO,  A  H°mo  =  227.5  kJ/mol  (3) 
CO  +  H20  =  H2  +  C02 ,  A  H°li00  =  -31.8  kJ/mol  (4) 

The  reaction  given  by  Eq.  (3)  is  a  slow  and  highly 
endothermic  reaction,  and  the  reaction  given  by  Eq.  (4)  is 
a  fast  and  weak  exothermic  reaction  [5].  The  latter  reaction, 
the  water  gas  shift  reaction,  is  assumed  to  be  in  equilibrium 
at  the  reforming  temperature. 


4.1.1.  Pre-reforming  of  methane 

In  pre-reforming,  carbon  formation  may  take  place  under 
unfavorable  conditions.  The  dominant  reactions  for  the 
carbon  formation  are  the  following  reactions  [6,7]: 

2CO  =  C02  +  C  (Boudouard),  (5) 

CO  +  H2  =  H20  +  C  (CO  hydrogenation) ,  (6) 

CH4  =  2H2  +  C  (CH4  decomposition) .  (7) 


The  activity  of  carbon,  Ac,  of  Eq.  (5)  is  given  a 

A ,  =  ^. 

pco2 


(8) 


The  activities  of  carbon  in  Eqs.  (6)  and  (7)  can  be  given  in 
the  same  manner.  These  activities  are  equal  in  equilibrium. 
In  nonequilibrium,  the  lowest  activity  has  the  priority  to  be 
checked.  Activity  exceeding  1.0  indicates  the  initiation  of 
carbon  formation.  As  mentioned  above,  water  vapor  neces¬ 
sary  for  the  reforming  and  the  suppression  of  carbon  for¬ 
mation  is  supplied  by  the  recirculation  of  the  spent  fuel.  The 
changes  in  the  carbon  activity,  Ac,  and  the  methane  reform¬ 
ing  rate,  Rr  (equal  to  the  ratio  reformed  CFLi/input  CFI4), 


with  the  reforming  temperature  and  the  fuel  recirculation 
rate,  f?ecf,  are  shown  in  Fig.  2.  Both  the  reforming  rate  and 
the  carbon  activity  are  sensitive  to  the  temperature  and  the 
fuel  recirculation  rate. 

The  pre-reforming  temperature  of  the  developed  SOFC 
was  973  K  [8].  As  can  be  seen  in  Fig.  2,  at  the  fuel 
recirculation  rate  of  0.4,  carbon  formation  may  take  place 
at  970  K.  This  temperature  has  no  margin  for  the  above- 
mentioned  pre-reforming  temperature.  Therefore,  the  fuel 
recirculation  rate  is  set  at  0.5  at  which  the  carbon  formation 
temperature  is  900  K. 

The  pre-reforming  rate  of  CH4  has  to  be  smaller  than  the 
value  given  in  Fig.  2  due  to  the  slow  reforming  reaction  and 
the  limit  of  the  pre-reformer  size. 


4.1.2.  Internal  reforming 

After  appropriate  pre-reforming,  the  fuel  containing  resi¬ 
dual  methane  goes  into  the  internal  reformer.  The  rate  of  the 
reforming  reaction,  f?CH4,  of  Eq.  (3)  is  given  as  [9] 


Rciu 


1.75  exp 


-57840\ 
RT  ) 


(9) 


The  change  in  the  methane  molar  flow,  dMCH4/dr,  is 
given  as 


dA7(  ■[  ]. 
dr 


-ArefflaqWcat. 


(10) 


where  Aref  is  the  reforming  area  per  unit  length  of  the 
internal  reformer  (or  cell),  and  WCSLt  the  catalyst  weight 
density. 

With  the  advance  of  the  methane  reforming  as  given  by 
Eq.  (10),  H2  and  CO  are  generated  consuming  CH4  and  H20 
as  given  by  Eq.  (3).  At  the  same  time,  the  water  gas  shift 
reaction  given  by  Eq.  (4)  reaches  equilibrium  rapidly.  There¬ 
fore,  CO,  C02,  FI2  and  HzO  have  to  satisfy  the  following 
equilibrium  equation: 


XcqAh2o  ’ 


(ID 


where  Ks  is  the  equilibrium  constant  of  Eq.  (4).  The 
atomic  conservation  of  C,  H  and  O  has  also  to  be  kept 


among  C02,  H2  and  H20,  excluding  CH4,  which  is  not  in 
equilibrium. 

4.2.  Cell 

4.2.1.  Equivalent  electrical  circuit  of  cell 

The  segmented  type  tubular  cell  stack  consists  of  many 
pairs  of  cells  and  interconnectors.  Fig.  3  shows  a  part  of  the 
cross-section  of  the  stack.  The  corresponding  equivalent 
ladder  circuit  of  the  cell  area  is  shown  in  Fig.  4. 

In  the  connecting  area,  a  similar  circuit  without  the  local 
electromotive  force  can  be  drawn  to  calculate  the  voltage 
drop.  Each  resistance  and  the  related  values  of  typical  cell 
components  at  1273  K  are  given  in  Table  1.  These  values 
were  obtained  from  the  literature  [10]. 

The  resistivity  of  the  electrolyte  strongly  depends  on  the 
temperature,  and  its  effect  cannot  be  ignored.  The  depen¬ 
dence  of  resistivity  on  the  temperature  is  given  by  the 
following  equation  [11]: 

Re  =Re0  exp  [lOlOO (12) 

The  resistivities  of  other  parts  are  assumed  to  be  constant 
because  these  resistivities  depend  weakly  on  temperature 
[12],  and  their  contributions  to  the  total  voltage  drop  are  small. 

Each  branch  current  can  be  obtained  by  the  iterative 
method,  giving  an  electromotive  force  E  and  a  polarization 
voltage  Vp.  The  method  of  calculating  the  electromotive 
force  and  the  polarization  will  be  given  in  the  following 
sections. 

Finally,  the  cell  voltage  and  then  the  output  voltage  of  the 
stack  can  be  obtained.  This  enables  calculating  the  effi¬ 
ciency  by  Eq.  (1). 


Fig.  4.  Equivalent  ladder  circuit  of  cell  area. 


4.2.2.  Changes  in  gaseous  component,  electromotive  force 
and  concentration  polarization  in  cell  stack 

The  changes  in  molar  fractions  of  the  fuel  in  the  cell  stack 
have  to  be  calculated  considering  the  reforming  of  residual 
methane  and  the  inflow  of  oxygen  due  to  the  cell  reaction. 

In  this  region,  CO,  C02,  H2  and  H20  have  to  satisfy  the 
equilibrium  equation  given  in  Eq.  (11).  However,  in  atomic 
conservation,  an  increase  in  O  due  to  the  cell  reaction  has  to 
be  considered. 

The  ideal  electromotive  force  of  the  relevant  position,  E, 
is  given  as 


E  =  E0(T)+—  In 


Xh2o 


(13) 


The  concentration  polarization  is  the  difference  between 
the  ideal  local  electromotive  force  and  the  real  electromotive 
force.  The  real  local  electromotive  force  can  be  given  in 
the  same  manner  as  Eq.  (13)  using  molar  fractions  at  the 
interfaces  between  the  electrode  and  the  electrolyte.  The 
concentration  polarization  of  the  relevant  position,  Vcons,  is 
then  given  as 


RT  2fH2ArH2o,B  J  X'Q 

Vcons  =  In - V - ,  (14) 

2F  XHlfiXH,o^X>0iR 

These  molar  fractions  also  have  to  satisfy  Eq.  (11).  To 
obtain  Vcons,  the  mass  fluxes  in  the  electrodes  and  the 
substrate  have  to  be  calculated. 

In  porous  media,  only  radial  flow  is  assumed  at  the  active 
cell  area,  and  no  flow  is  assumed  in  the  interconnecting  area. 


Table  1 

Resistance  and  diffusion  parameters  at  1273  K  [10]a 


Component 

Thickness  (cm) 

Resistivity  (12  cm) 

Resistance  (12) 

Air  electrode  (Sr  doped  LaMn03) 

0.1 

0.013 

Ra  =  0.0184(5* 

Electrolyte  (Y203  doped  Zr02:  YSZ) 

0.005 

10.0 

Re  0  =  0.00745/5* 

Fuel  electrode  (40  vol.%  Ni/YSZ  cermet) 

0.015 

0.001 

Rf  =  0.010035* 

Interconnector  (Mg  doped  LaCr03) 

0.003 

0.5 

Rint  =  0.00022/5* 

a  Total  stack  length  =108  cm;  unit  cell  length  =  1.2  cm  active  area  plus  0.6  cm  connecting  area;  substrate:  outer  diameter  D  =  2.1  cm,  thickness  = 
0.2  cm;  porosity  e  =  0.4  (common  in  electrodes  and  substrate). 
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Here,  we  define  the  current  density  J  as  the  current  density  at 
the  surface  of  the  substrate 


2nrsvlbdx ' 


(15) 


The  molar  fluxes  of  the  gas  components  in  porous  media 


have  the  following  relations: 

Vh,+«co=(^)V, 

(16a) 

Nh2o  =  — Nh2,  Nco2  =  —Ncch 

(16b) 

(16c) 

The  effect  of  the  direct  oxidation  of  methane  to  contribute 
to  Eq.  (16)  is  neglected  because  its  oxidation  speed  is  slow. 
The  gaseous  flows  in  the  porous  media  are  assumed  to  be 
governed  by  the  Stefan-Maxwell  equations  as  [13] 

iLY,  XiNj-XjNj 

dr  cDy  ’ 

(17) 

Dtj  is  given  as  [14] 

Djj  =  Dp2, 

(18) 

where  Dps  the  diffusivity  between  gas  species  i  and  j  in  free 
space  that  can  be  calculated  by  the  Chapman-Enskog 
equation  [15]  and  s  the  porosity  of  electrode  or  substrate. 

The  total  molar  density,  c,  is  assumed  to  be  constant  at  any 
position.  The  molar  fractions  at  the  main  flow  and  the 
interface  can  be  calculated  using  Eqs.  (9)— (1 1),  (15)— (18). 
We  can  then  calculate  the  concentration  polarization  by 
Eq.  (14). 


4.2.3.  Activation  polarizations 

The  activation  polarization,  Vact,  was  given  by  Sawata 
et  al.  as  follows:  [16] 

where 


_  oqRT 

0  “  (2  + a )F 


(20) 


where  a  is  the  apparent  transfer  coefficient  denoting  unsym- 
metrical  reaction,  a0  the  interface  conductivity.  Setoguti 
et  al.  [17]  suggested  that  a0  for  the  fuel  electrode  mainly 
depends  on  the  residual  oxygen  pressure  in  the  fuel.  How¬ 
ever,  that  relation  is  complicated.  Here  we  simplify  the 
dependence  of  a0  on  the  oxygen  pressure,  Po2,b,  as  follows: 


<io  =  P  exp  | p'o2  .h  (21) 

where  /]  is  a  constant,  y  the  reaction  order,  AGact  the 
activation  energy.  In  the  oxidant,  Pq2,b  =  P'q  b-  In  the  fuel. 


Table  2 

Coefficients  to  give  activation  poianzauon 

Position  a  />’  y  AGact  (kj/mol) 

Fuel  electrode  1  107'79  0.133  120 

Air  electrode  2  10602  0.5  130 


the  oxygen  pressure  is  given  as 


(22) 


Yamamoto  et  al.  [18]  suggested  that  y  for  the  air  electrode 
was  approximately  0.5. 

From  our  small  experiments  and  these  references,  each 
coefficient  is  given  as  shown  in  Table  2. 

The  activation  polarization  Vactox  in  the  oxidant  and  the 
activation  polarization  Vact  f  in  the  fuel  can  be  calculated 
numerically  using  Eqs.  (19)-(22).  The  total  polarization  is 
then  given  as 


Vp  =  Vcons  +  Vact, ox  +  Vactf. 


(23) 


4.2.4.  Total  voltage  drop  and  cell  voltage 

Now  we  can  calculate  the  current  distribution  and  then  the 
voltage  drop  in  the  cell  stack.  The  average  voltage  drop  and 
the  average  cell  voltage  are  shown  in  Fig.  5.  The  simulated 
cell  voltage  of  this  figure  agrees  well  with  that  reported  by 
Singhal  [19]  as  seen  in  Fig.  5. 

4.3.  Enthalpy  balance 

The  enthalpy  balance  is  calculated  using  the  control- 
volume  method  [20],  The  term  control-volume  is  abbreviated 
as  CV  here.  The  dimension  of  the  specified  control-volume 
along  flow  direction  is  taken  as  3x. 


1273  K.  Eav:  average  electromotive  force;  Fceu:  average  cell  voltage;  Vine: 
average  ohmic  voltage  drop  in  interconnector;  V7uei:  average  ohmic  voltage 
drop  in  fuel  electrode;  Vi,,,:  average  ohmic  voltage  drop  in  air  electrode; 
Vdecr  average  ohmic  voltage  drop  in  electrolyte;  WH:  V^n  reported  by 
Singhal  of  SWHPC. 
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Fig.  6.  Enthalpy  balance  in  CVs  of  solid  and  gas  phases. 


All  solids  are  assumed  to  be  uniform  solids  having  the 
equivalent  thermal  conductivities  of  multilayers. 

All  the  chemical  reactions  are  assumed  to  take  place  in  the 
solids  (catalyst,  electrodes,  interfaces  of  electrode  and  elec¬ 
trolyte).  The  schematic  heat  migration  in  the  solid  CV  and 
the  gas  CV  is  shown  in  Fig.  6.  Each  CV  is  numbered 
sequentially  as  C V(i  -  1),  CV(i),  CV(i  +  1).  Here,  the  heat 
migration  from  CV(i  —  1)  to  CV(i)  due  to  heat  conduction  is 
denoted  as  GcondO  —  1),  the  heat  migration  due  to  the 
convection  is  denoted  as  QcnnJi  —  1),  and  so  on. 

4.3.1.  Enthalpy  balance  in  gas  CV 

The  increase  in  the  enthalpy  in  CV  lying  between  x  and 
x  +  bx  has  to  balance  the  enthalpy  input  from  the  side  walls. 
For  the  gas  CV(i)  (see  Fig.  6), 


Gconv(i)  -  Qconv(l  ~  1)  =  GLd  +  E  ^ans  "  Gre 


(24) 

where 

Gconv(i)  =  EM7ff7(r g(0), 

(25) 

rPg(i) 

Hj(T$(i))=H?(T0)+  CPjdT. 

Jto 

(26) 

Each  enthalpy  on  the  right  side  of  Eq.  (24)  is  given 
follows: 


QL d  =  ln(‘  +  1)  +  Tg(i  -  1)  -  2Tg(0]  ^ ,  (27) 

Gtnms  =  -  Tg(i%  (28) 

Greact  =  Greform  +  Gshift, 

(for  the  fuel  in  the  internal  reforming  region)  (29a) 
Greact  =  Gceii.ox ,  (for  the  oxidant  in  the  cell  region) 

(29b) 

Greact  =  Greform  +  Gshift  +  Gccll.f, 

(for  the  fuel  in  the  cell  region)  (29c) 

Greact  =  0,  (for  other  region)  (29d) 


where 

Greform  =  SMCHi[HCHi(Tg(i))  +  HHl0(Ts(i))  -  3//Hj  (Ts(l)) 
-Hco(Tsm,  (30) 

where  c5/Wch4  moles  of  reformed  CH4  in  solid  CV(i)  that  can 
be  calculated  using  Eq.  (10), 

Gshift  =  SMC02[HC0(Ts(i))  +HH20(Ts(i))-Hn2(Ts(i)) 

~  HC02(Ts(t)%  (31) 

where  5MCo2  is  the  moles  of  C02  generated  by  the  water  gas 
shift  reaction  in  solid  CV(i)  that  can  be  calculated  using 
Eq.  (11)  and  atomic  conservation.  If  5Mco2  is  negative,  T„(i) 
and  Ts(i)  in  the  second  term  on  the  right-hand  side  of  Eq.  (3 1 ) 
have  to  be  exchanged: 

Gcell,ox=^ffo2(TOx(0),  (32) 

Gccii.f  =  dMn2[HHl(Tf(i))  -  //1,2o(rccii(i))J 

+  SMCo[Hco{Tf(i))  -  HCo2(TceU,f(i))],  (33) 

where  dM^2  and  SMqo  can  be  obtained  using  Nu2  and  Nqo  as 
solutions  for  Eqs.  (16)  and  (17). 

4.3.2.  Enthalpy  balance  in  solid  CV 

The  enthalpy  balance  in  CV  of  solid  is  given  as 

Gcond(0-Gcond(;  -  1)  =  E  Grad  +  Greact^ ~E  GOars  “ Ggene , 

(34) 

where 


2cond(0  —  ^sks 


rg(i+i)-7’,(Q 

Sx 


(35) 


and  Grad  is  the  enthalpy  input  due  to  radiation  from  neigh¬ 
boring  solid  CVs  [21],  Greact  and  G^ans  the  same  as  those 
defined  in  the  preceding  section  and  Ggene  the  cooling  due  to 
power  generation  for  the  cell. 

Each  enthalpy  on  the  right-hand  side  of  Eq.  (34)  is  given 
as  follows: 


Grad 

5.67(7^ -7f(i))10-12 
"  l/(es5sl)  +  (1/ (es  -  l))/5s2  ’ 

(36) 

Greae 

t  =  Greform  +  Gshift,  (for  the  internal  reformer) 

(37a) 

Greae 

t  =  Greform  +  Gshift  +  Gcell.f  +  Gcell.ox, 

,  (for  the  cell) 

(37b) 

Greae 

t  =  0,  (for  the  solid  parts) 

(37c) 

Ggene 

,  =  0,  (except  for  the  cell  part) 

(38a) 

Ggene  =  4(*)[ff(0  ~  Vp (i)  -  /e(t)*e(0] 

[ll{l)R&  +  42(*  -  l)/?a  +  /f(i)/?f  +  If(i  -  l)/?f] 

2 

(for  the  cell)  (38b) 
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In  the  connecting  area,  Vp(i )  vanishes  in  Eq.  (38a),  and  the 
currents  and  resistances  have  to  be  replaced  by  appropriate 
values  matching  the  position. 

Since  the  Reynolds  numbers  of  die  fuel  and  the  oxidant 
are  under  500,  die  heat  transfer  coefficient,  hg,  is  calculated 
using  a  fully  developed  laminar  flow  approximation  at 
constant  wall  temperature  as  [22] 


hg 


3.66kg 

de 


(39) 


where  kg  is  the  thermal  conductivity  of  the  mixed  gas,  de  the 
equivalent  diameter,  kg  the  calculated  using  the  Bromley- 
Wike  equation  [23].  The  thermal  conductivities  of  gases  are 
taken  from  the  data  book  [24],  The  thermal  conductivities  of 
the  solid  components  are  given  in  Table  3.  The  dependence 
of  these  values  on  temperature  is  ignored  for  simplicity. 

The  calculated  typical  temperature  distributions  are 
shown  in  Fig.  7.  The  temperature  distribution  of  the  cell 
stack  is  moderate,  and  the  average  temperature  is  1260  K. 
There  is  no  literature  to  refer  to  concerning  this  temperature 
distribution.  However,  the  tendency  of  the  temperature 
distribution  of  the  cell  stack  agrees  well  with  that  of 
3  kW  SOFC  of  WH  [27], 


4.3.3.  Exhaust  temperature 

The  exhaust  temperature  can  be  obtained  from  the  follow¬ 
ing  global  heat  balance: 

2  MjHj(Tmt)=  ]T  M'-Hj{TAQ)  +  AfCH4tfcH4(7>o)  -Vo 

(40) 

For  the  input  air, 

,  Ur  ,  ( 0.79\  , 

M'0i  =  2 MCHi  yL,  M'N2=  ^0  2|J M'Qi ,  (41) 

where  Ar,  C02  and  other  elements  are  neglected  as  compo¬ 
nents  of  air. 

For  the  exhaust  gas, 

Mq2  =  2MCh4  (Jpj  I-  MNl=  M'N2  (42) 

Mco2  =  McUi,  Mh2o  =  2MCh4-  (43) 


5.  Procedure  for  numerical  results 

The  main  results  desired  are  the  current  distribution,  the 
gas  component  distribution  and  the  temperature  distribution. 
Each  computation  is  carried  out  using  die  corresponding 
models.  An  algorithm  to  obtain  numerical  results  is  given  in 
Fig.  8.  Each  computation  step  is  executed  using  the  results 
given  in  the  former  computation  steps.  The  cycle  of  the 
computation  steps  is  repeated  until  global  convergence  is 
obtained. 

The  effect  of  CV  size  on  the  computation  results  is  given 
in  Table  4.  However,  the  size  has  little  effect  on  the  results. 
Therefore,  the  CV  size  is  taken  at  0.3  cm  for  the  speed  of 
convergence  of  the  computation. 


6.  Specification  of  operating  conditions  of 
SOFC  system 

The  fuel  utilization,  Uf,  is  set  at  85%  referring  to  the 
reported  experiments  as  mentioned  in  Section  3. 

The  target  value  of  the  efficiency  of  an  SOFC  operating  at 
1273  K  is  about  50%  (HHV)  [28],  Therefore,  Vcell  has  to  be 
about  0.68  V  as  deduced  from  Eq.  (1).  As  seen  in  Fig.  5,  the 
cell  voltage  of  0.68  V  can  be  obtained  at  a  current  density  of 
about  0.3  A/cm2;  therefore,  the  current  density  is  set  at 
0.3  A/cm2. 

The  inlet  temperature  of  the  raw  fuel,  Tm  (see  Fig.  1), 
is  preset  at  288  K.  The  inlet  temperature  of  the  air,  T.d0 
(see  Fig.  1),  is  the  same.  This  temperature,  288  K,  is  the 
representative  ambient  temperature  adopted  in  gas  turbine 
systems  [29], 

The  pre-reforming  temperature  (equal  to  preheated  tem¬ 
perature  of  the  fuel),  rfln  (see  Fig.  1),  is  set  at  953  K.  This 


Fig.  8.  Algorithm  to  obtain  numerical  results. 


temperature  is  the  value  which  is  approximate  to  the  reported 
reforming  temperature  [8]  and  avoids  the  formation  of 
carbon.  The  preheated  air  temperature,  (see  Fig.  1),  is  set 
at  993  K.  By  setting  these  preheated  temperatures,  the 
temperature  difference  between  the  maximum  temperature 
and  minimum  temperature  can  be  limited  to  within  1 13  K  as 
seen  in  Fig.  7. 

The  operating  pressure  of  the  cell  stack  is  set  at  1.0  atm. 
Though  it  is  necessary  to  carry  out  a  slight  pressurization, 
this  effect  is  neglected  for  simplicity. 

Table  4 

Effect  of  CV  size  on  results 

(cm)  %c  (%)  Um  (%)  Tout  (K) 

48.397  49.961  959.45 

5  48.455  49.938  958.55 


The  fuel  recirculation  rate  is  set  at  50%  to  suppress  carbon 
formation  as  described  in  Section  4.1.1.  The  oxidant  recir¬ 
culation  rate  is  also  set  at  50%,  but  there  is  no  reason  to  set 
this  value.  However,  the  effect  of  recirculation  will  be 
discussed  later. 

The  internal  reforming  rate  is  set  at  90%  (equal  to  10% 
pre-reforming)  in  order  to  increase  the  cooling  effect  by  the 
internal  reforming.  For  this  reforming  rate,  the  catalyst 
density  is  set  to  2.3  mg/cm2  by  trial  and  error.  The  pre¬ 
reforming  of  10%  is  a  value  considering  the  progress  of 
some  reforming  in  the  mixing  of  raw  fuel  and  recycled  fuel. 

These  operating  conditions  are  chosen  to  obtain  a  mod¬ 
erate  temperature  distribution  in  the  cell  stack.  The  effects  of 
these  operation  conditions  of  an  SOFC  system  will  be 
discussed  later. 

The  maximum  temperature  of  the  cell  stack  is  set  at 
1300  K  by  adjusting  the  oxygen  utilization  rate. 

With  this  maximum  temperature  and  these  system  oper¬ 
ating  conditions,  we  can  obtain  the  average  cell  temperature 
of  1260  K  which  is  almost  equal  to  the  1273  K  for  the  state- 
of-the-art  of  the  operation  temperature  of  the  SOFC  [27]. 


7.  Results  of  simulation  and  discussion 

7.1.  Effect  of  catalyst  density 

By  setting  the  appropriate  uniform  catalyst  density  in  the 
internal  reformer,  the  reforming  rate  reaches  approximately 
95%  at  the  outlet  of  the  internal  reformer,  as  shown  by  curve 
(B)  in  Fig.  9.  In  the  case,  the  temperature  distribution  of  the 
cell  is  moderate,  as  shown  by  curve  (b)  in  Fig.  10. 

In  the  case  of  a  thin  catalyst  density,  the  reforming  rate  at 
the  outlet  of  the  internal  reformer  is  not  sufficient  as  shown 
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Fig.  10.  Effect  of  catalyst  density  on  cell  temperal 


by  curve  (A)  in  Fig.  9.  Therefore,  much  methane  not  yet 
reformed  remains.  This  methane  is  reformed  quickly  at  the 
fuel  inlet  region  (at  the  right  end  of  the  cell)  decreasing  the 
cell  temperature  as  shown  by  curve  (a)  in  Fig.  10,  because 
the  fuel  electrode  acts  as  a  thick  catalyst  for  reforming. 

On  the  contrary,  in  the  case  of  a  thick  catalyst  density,  the 
reforming  occurs  quickly  at  the  inlet  region  of  the  internal 
reformer,  decreasing  the  neighboring  cell  temperature  as 
shown  by  curve  (c)  in  Fig.  10.  This  decrease  in  the  tem¬ 
perature  increases  the  carbon  activity,  Ac,  to  nearly  1.0  as 
shown  by  curve  (c)  in  Fig.  9. 

The  catalyst  density  also  affects  the  efficiency,  the  oxygen 
utilization  and  the  exhaust  temperature.  The  maximum 
exhaust  temperature,  Tout,  and  oxygen  utilization,  Uox,  can 
be  obtained  at  a  moderate  catalyst  density  as  shown  in 
Fig.  11.  However,  the  effect  on  the  efficiency  is  insignificant. 

7.2.  Effect  of  internal  reforming 

In  this  section,  the  effect  of  the  internal  reforming  rate  is 
discussed  with  the  relative  catalyst  density  one  given  in  the 


Fig.  1 1 .  Effect  of  catalyst  on  output  characterist 


preceding  section.  As  mentioned  above,  the  reforming  is  a 
strongly  endothermic  reaction.  Therefore,  the  increase  in  the 
internal  reforming  rate  needs  a  decrease  in  the  air  flow  rate 
(namely,  an  increase  in  the  oxygen  utilization  t/ox)  to  keep 
the  total  cooling  capacity  at  the  same  level.  As  a  result,  the 
increase  in  the  internal  reforming  rate  decreases  the  electro¬ 
motive  force  and  hence  decreases  the  efficiency,  as  shown  in 
Fig.  12.  The  decrease  in  the  air  flow  increases  the  exhaust 
temperature. 

7.3.  Effect  of  graded  catalyst  distribution 

In  the  previous  sections,  the  catalyst  density  was  uniform 
throughout  the  internal  reformer.  In  this  section,  the  effect  of 
the  graded  catalyst  density  will  be  discussed.  The  rate  of 
internal  reforming  is  kept  constant  at  90%. 

The  total  area  of  the  internal  reformer  is  divided  into 
several  regions  of  equal  length.  The  setting  of  the  catalyst 
density  for  each  region  is  so  calculated  as  to  obtain  the 
maximum  efficiency,  while  the  total  catalyst  weight  is  kept 
constant  and  the  gradient  of  the  density  distribution  is  the 
upper  right  ascent.  The  obtained  catalyst  density  for  each 
region  depends  on  the  number  of  divisions  as  shown  in 
Fig.  13. 

The  effect  of  the  number  of  the  regions  on  the  cell 
temperature  distribution  is  shown  in  Fig.  14.  The  curve 
for  one  region  is  omitted  because  it  is  the  same  as 
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Fig.  13.  Distribution  of  graded  catalyst  density.  Internal  reforming,  90%. 
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Fig.  10(b).  As  shown  in  the  figure,  the  increase  in  the  average 
temperature  of  the  cell  stack  with  the  increasing  number  of 
regions  can  be  seen.  This  increase  in  the  average  temperature 
decreases  the  internal  losses  and  therefore  decreases  the  air 
flow  rate  for  cooling  the  cell.  However,  this  decrease  in  the 
air  flow  rate  decreases  the  electromotive  force.  As  a  result, 
the  efficiency  and  the  exhaust  temperature  slightly  increase 
as  shown  in  Fig.  15. 

7.4.  Effect  of  air  inlet  temperature  and 
oxidant  recirculation 

Here,  the  effect  of  each  operating  condition  specified  in 
Section  6  on  the  output  characteristics  will  be  examined.  As 
mentioned  in  Section  2,  part  of  exhaust  oxidant  is  recircu¬ 
lated  to  heat  the  inlet  air  to  Tain.  The  effect  of  the  air  inlet 
temperature,  Tain,  is  shown  in  Fig.  16. 

The  air  (oxidant)  functions  as  a  coolant.  Therefore,  the 
increase  in  the  air  inlet  temperature  needs  an  increase  in  the 
air  flow  rate  (namely,  a  decrease  in  the  oxygen  utilization 
Uox)  to  keep  the  cooling  capacity  at  the  same  level.  The 
increased  air  flow  dilutes  the  combustion  exhaust  gas  to 
decrease  the  exhaust  temperature.  On  the  other  hand,  the 
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Fig.  16.  Effect  of  air  inlet  temperature. 
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increase  in  the  air  flow  rate  increases  the  output  voltage  and 
hence  increases  the  efficiency.  However,  this  effect  is  insig¬ 
nificant. 

The  effect  of  the  oxidant  recirculation  is  shown  in  Fig.  17. 
Here,  T is  kept  constant  at  993  K.  The  increase  in  the 
oxidant  recirculation  rate  needs  an  increase  in  the  oxidant 
utilization  to  keep  the  cooling  capacity  at  the  same  level.  The 
increase  in  the  oxidant  recirculation  leads  to  a  decrease  in 
the  oxygen  partial  pressure  in  the  cell  stack.  Therefore,  the 
increase  in  the  oxidant  recirculation  decreases  the  cell. 

7.5.  Effect  of  raw  fuel  inlet  temperature  and 
fuel  recirculation 

The  effect  of  the  inlet  temperature,  Tfin,  of  the  partially 
pre-reformed  fuel  is  shown  in  Fig.  18.  The  increase  in  the 
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Fig.  19.  Effect  of  fuel  recirculation. 

inlet  temperature  decreases  the  cooling  capacity.  Therefore, 
the  flow  rate  of  the  air  increases  to  compensate  for  the 
cooling  capacity.  The  effect  of  the  increase  in  the  fuel  inlet 
temperature  resembles  that  of  the  increasing  air  inlet  tem¬ 
perature.  However,  it  is  less  effective  because  the  heat 
capacity  of  the  fuel  is  smaller. 

The  effect  of  the  fuel  recirculation  is  shown  in  Fig.  19.  In 
this  case,  Tfin  is  kept  constant  at  953  K.  To  avoid  carbon 
formation,  water  is  added  at  the  lower  fuel  recirculation. 
Therefore,  the  total  cooling  capacity  of  the  fuel  and  the 
recirculated  fuel  is  kept  almost  constant.  Hence,  the  oxidant 
utilization  only  changes  slightly.  The  efficiency  hardly 
changes  with  the  fuel  recirculation  rate,  because  the  oxidant 
flow  rate  and  the  H2/H20  ratio  in  the  fuel  are  kept  at  almost 
constant  (namely,  the  electromotive  force  hardly  changes). 
The  heat  necessary  for  the  vaporization  of  the  added  water  is 
supplied  from  the  combustion  exhaust  gas  in  the  pre-refor¬ 
mer.  Therefore,  the  exhaust  temperature  from  the  SOFC 
system  decreases  with  the  decrease  in  the  fuel  recirculation. 

8.  Summary  and  conclusion 

The  effect  of  the  internal  reforming  on  the  output  char¬ 
acteristics  of  the  SOFC  system  was  discussed  using  numer¬ 
ical  models.  Some  effects  of  gas  inlet  conditions  have  also 
been  discussed.  The  following  conclusions  have  been  drawn 
for  several  sets  of  operating  conditions  of  the  SOFC  system 
without  the  heat  leak,  restricting  the  maximum  cell  tem¬ 
perature  at  a  constant  1300  K  by  adjusting  the  oxidant 
utilization  rate. 

1.  Setting  an  appropriate  uniform  catalyst  density  leads  to 
a  moderate  temperature  distribution  and  leads  to  a 
higher  temperature  of  the  exhaust  gas. 

2.  Increasing  the  internal  reforming  rate  leads  to  a  large 
increase  in  the  exhaust  temperature  and  leads  to  a  slight 
decrease  in  the  efficiency. 

3.  Setting  the  graded  catalyst  density  leads  to  a  flatter 
temperature  distribution  in  the  cell  stack  and  leads  to  a 
slight  increase  in  the  exhaust  temperature. 


4.  Increasing  the  inlet  gas  temperature,  leads  to  a  slight 
increase  in  the  efficiency;  however,  leads  to  a  decrease 
in  the  exhaust  gas  temperature  from  the  SOFC  system. 
The  effect  of  the  air  inlet  temperature  on  the  exhaust 
temperature  is  significantly  stronger  than  that  of  the  fuel 
inlet  temperature. 

5.  Increasing  the  oxidant  recirculation  leads  to  a  slight 
decrease  in  the  efficiency  but  leads  to  a  large  increase  in 
the  exhaust  gas  temperature.  As  for  the  fuel  recircula¬ 
tion,  the  tendency  is  similar. 

These  results  can  be  summarized  as  the  operating  con¬ 
ditions  of  the  present  SOFC  system  which  have  considerable 
effects  on  the  exhaust  temperature  but  a  slight  effect  on  the 
efficiency. 

The  effects  of  the  heat  leak  and  the  pressurization  will  be 
investigated  in  the  future  in  the  relation  with  the  combined 
power  generation  cycle. 
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